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Abstract: We study the infrared (large separations) behavior of a massless minimany 
coupled scalar QFT with a quartic self-interaction in de Sitter. We show that the perturba- 
tion series in the interaction strength A is secular, i.e it does not lead to a uniform approx- 
imation of the solution in the infrared region. Only a non-perturbative resummation can 
capture the correct IR behavior. We seek to justify this picture using the Dyson-Schwinger 
equations in the ladder-rainbow approximation. We are able to write down an ordinary 
differential equation obeyed by the two-point function and perform its asymptotic analysis. 
Opposing a widespread belief that IR divergences are regulated through a dynamical mass 
generation, our work suggest an even stronger IR decay than the power-law one resulting 
from a mass term. 
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1 Introduction 



1.1 Prelude 

The two-point function A^2 of a massive scalar field in de Sitter decays in the IR region, 
i.e for largely separated spacetime points. For small masses, A^2 develops a 1/m? pole, 
making the massless limit ill-defined and the quantization of the theory a non-trivial task. 
This situation is to contrast with the flat space case, where massive and massless fields two- 
point functions decay for large invariant distance /i (as //~^e~"^^^ and ;U~^ respectively). 

Several propositions to quantize the massless minimally coupled (mmc) field co-exist 
in the literature, depending on how one interprets and cures this ill-defined massless limit. 
The first and most popular direction is to understand the problematic massless limit as an 
indication of the impossibility of maintaining the full de Sitter invariance of the quantum 
theory. Indeed, breaking de Sitter invariance leads to a well-defined quantization and to 
interesting cosmological implications [1, 2]. 

Another option, investigated in [3, 4] (see also the earlier Starobinsky's stochastic 
approach [5, 6]), is to recognize the necessity of the inclusion of a small (say quartic) 
interaction of strength A. The idea being that the presence of a small but non-vanishing 
A will make the massless limit smooth (in other terms, the limits m — t- and A — ?• 
do not commute). These investigations lead to the conclusion that interactions induce a 
dynamically generated mass m^g oc \/Aif-where H is the Hubble constant-that regulates 
the IR behavior. In particular, in [4] an interesting exact treatment of the zero-mode in 
the interacting theory (which amounts to performing a non-gaussian but one dimensional 
integral) is given, and the author argues in favor of a constant massless two-point function 
proportional to This approach, as recently shown in [7] leads again to the conclusion 

that interactions induce a dynamically generated mass. 

Yet another approach can be found in [8, 9], where a "renormalized" two-point function 
Aq is defined through the subtraction of the divergent 1 /vn? term. This approach leads to 
a well-defined two-point function Aq, which is growing in the IR. The construction of Aq 
presents a lot of advantages, as it preserves causality and positivity. We have also shown in 
[10] that it leads to the observed scale-invariant CMB power-spectrum. Our present work 
is based on this free propagator. 

While we indeed believe that the inclusion of an interaction, whatever small, has 
drastic consequences on the quantization of the massless field and its IR behavior, we 
argue here in favor of a quite different overall picture. The free propagator Aq, as well 
as its perturbative corrections are growing. However, when the leading IR contributions 
are summed up, the full non-perturbative propagator G is decaying. In other words, the 
perturbative expansion develop secular terms, making the perturbative expansion non- 
uniform in the IR. Contrary to the previously discussed papers, we suggest that the decay 
of G is too rapid to be interpreted as resulting from a dynamical mass generation. We take 
a first step to justify this picture by summing, via a linear Dyson-Schwinger equation, the 
so called ladder-rainbow diagrams. 

The organization of the paper is as follows. After a rapid introduction to the basics of 
de Sitter geometry and QFT, we write down in section II the Dyson-Schwinger equations in 
the so-called ladder-rainbow approximation. Taking full-advantage of de Sitter symmetry. 
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we are able to reduce this integral equation into an ordinary differential equation. The 
asymptotic analysis of which is studied in section III. A concise discussion of dynamical 
mass generation is given in section IV. Finally, section V is a discussion of several important 
points. 

1.2 De Sitter geometry and QFT 

The d-dimensional de Sitter spacetime can be identified with the real one-sheeted hyper- 
boloid in the d+ 1 Minkowski spacetime M^_|_i: 

with R = H^^ > being the de Sitter "radius". This definition of the de Sitter manifold 
reveals the maximal symmetry of Xd under the action of the de Sitter group SOq{1, d). Let 
fi{x,x') denote the geodesic distance between the points x and x' . It is useful to introduce 
the quantity z given for spacelike separations {fi^ > or < z < 1) by z = cos^ (2^)- 
As a general rule, we will work in the spacelike region and then analytically continue the 
propagators. As shown in [11], when acting on functions of the invariant the 
Laplace-Beltrami operator reduces to the ordinary differential operator : 



— + {d-l)A{fi)- 
df-i^ djj. 

where 



-n = — + {d-l)A{^l)— (1.1) 




= , on Minkowski space 

= ^^^^ — \ on de Sitter space. 

In terms of the variable z, this becomes the hypergeometric operator 

d 

-D = H = z(l - z)—r + [c - (a + b + l)z] — - ah 
dz'^ dz 

where 



_d-l + v^(cf - 1)2 - 4m2i?2 d - 1 - ^(d - 1)2 - 4m2i?2 

a — ^ , — ^ , c — a/ 1. 

In particular, in the massless case we have 
1.3 The massless field 

The physical reason behind the appearance of strong IR effects in de Sitter can be simply 
understood: the rapid expansion of the spacetime dilate correlation patterns. After all this 
is the exact reason why a de Sitter inflationary phase in the early universe solves many 
problems of the hot big-bang model. These IR effects are naturally stronger for massless 
(and non-conformally invariant) fields such as the massless minimally coupled (mmc) ^ and 
the graviton. We review here these IR divergences in the mmc case. 



^This means that the there is no term in the action proportional to TZ4P , TZ being the Ricci scalar, which 
is constant in de Sitter spacetime 
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Recall that in the Bunch-Davies vacuum state the massive propagator for a scalar field 
reads : 



A(x,x') =r -0-) \- V 2F1 



-a,a + d- l,-;z 



(1.3) 



where f = ~f + Y^~ m?E? and 2F1 is the hypergeometric function. In the massless 
limit, o" — 7- 0, this expression diverges because of the pole of the first Gamma function and 
we have the small mass expansion: 

A. r. d ^fd-l\ 1 

A(x, X ) ~ — STT^ ( — — ) — + regular terms m m. 
A-K' 



Note that in the flat space limit [R — )• 00), this singular term is absent and the massless 
limit is smooth 2. More precisely we have (for d> 2): 



where ^ = {x — x')^ is the invariant distance and i^^^) jg Hankel function of the 
second kind. 

One of the first papers studying the mmc scalar field in de Sitter is [1], where the 
authors prove that a usual de Sitter-invariant Fock space quantization is impossible in this 
case. They then propose to trade the de Sitter SO(l,(i) invariance for a smaller one, say 
a SO{d) invariance. Equivalently, it is a common belief among workers in the field that 
the scale-invariant power-spectrum leads necessarily to a breakdown of de Sitter invariance 
and that some physical quantities might thus become time-dependent. Several authors 
later proposed different treatments of the mmc field, among which [9] is one of the most 
exhaustive. Here the divergent term is subtracted, a "renormalized" de Sitter-invariant 
two-point function is computed ^ and it reads (we work in d = 4 for simplicity): 



Ao(^) ^ 



(47r)2i?2 



21n(l - z) 



(1.4) 



In [10] we also proved that Aq gives the observed scale-invariant power spectrum. We 
will use this propagator for the mmc hereafter (and we will omit the subscript when no 
confusion can arises). 



^This means that the flat space hmit (R — i> 00) and the massless limit (m 0) do not commute. This is 
a physically important fact and might mean that even a small amount of curvature - like in today's universe 
- might have important consequences on massless fields. 

^The draw-back is that the two-point function no longer verifies the equation of motion = 0, instead 
it verifies the anomalous equation: 

This simple renormalization procedure has been used implicitly in several earlier works. However, the 
major contribution of [9] is proving that on a suitably chosen subspace of states £, the equation of motion 
is effectively restored. This "Physical" space of states should be regarded the same way as we regard the 
one that appears in the quantization of gauge theory (for instance the space of transverse photons in QED). 
Moreover, the authors were able to show that the renormalized two-point function defines a positive kernel 
when restricted to £, thus enabling a probabilistic interpretation of the theory. 
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2 Dyson— Schwinger equations 



The Dyson-Schwinger equations, an infinite set of integral equations between the n-point 
functions, are the equations of motion of QFT. These equations, through some truncation 
schemes, give a convenient handle on non-perturbative effects in the theory. From now on 
we study the QFT given by the action: 



S = d x^J —g 



where (D/x = ^J—gd^x is the invariant volume element. 

The (unrenormalized) Dyson-Schwinger equations for the exact propagator involve the 
four-point function 

G{x,x') = A(x,x') + j dVy,dVy, A(x,yi)S(yi,y2) G{y2,x') (2.1) 
T:{x,y) = \{-iXf I dVz^dVz2dVz3G{x,zi)G{x,Z2)G{x,Z3)T4{zi,Z2,Z3,y) (2.2) 



where G is the exact Feynman propagator, A the free one, S the self-energy and r4 the 
exact 4-point function. A graphical representation of these two-equations is given in figure 
(1). The first equation is known as the Dyson equation which in flat spacetime is readily 
solved in Fourier space and is equivalent to the usual geometric series of self-energies: 



1 - A(p)S(p) 



6» 



Figure 1. Dyson-Schwinger equation relating the two and four point functions. 

Actually, let us now analyze the Dyson-Schwinger equations for our purposes in some 
detail. We concentrate on the flat space case, the results being immediately generalizable 
to de Sitter spacetime. To start with, we note that our self-energy starts at two loops. 
This is obvious, as also in de Sitter the two-point tadpole graph has zero renormalized 
value, the tadpole vanishes [12]. Thus, the expansion of the self-energy starts at two loops, 
with the first term in the expansion indeed furnishing a primitive element in the Hopf 
algebra H of cp'^ theory, which provides a Hochschild one-cocycle B^. Omitting higher 
Hochschild co-cycles -which seems reasonable as they do not alter the algebraic structure 
of the Dyson-Schwinger equations-, we find the combinatorial Dyson-Schwinger equation 
for the self-energy 

X(a2) = 1-^b/ ^ 



6 yx^x"^) 
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Figure 2. There is no contribution from this graph. 

Its fix-point is a formal series S i/[[A^]] such that the apphcation of renormalized Feynman 
rules gives ^r{X{X'^)) = 1 - S. 

Hochschild cohomology ensures that we renormalized by local counterterms, and en- 
sures that each contributing graph is divided by its symmetry factor, as it should [13, 14]. 

We drastically simplify this system by linearizing it to a commutative and co-commutative 
Hopf algebra generated from simple concatenations of the cocycle B^. 

X(A^) = 1 - (^(A^)) . 

Graphically, the two equations read as in Fig. 3. 
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Figure 3. Flat space non-linear and linear Dyson-Schwinger equations. 



Explicitly, the two expansions in graphs read as in Fig. 4. Note that in the lin- 



Figure 4. Non- linear and linear Dyson-Schwinger equations expanded out in perturbation theory. 
Note that the linear Dyson-Schwinger equations misses some graphs, which also results in wrong 
symmetry factors. It nevertheless gives a strict subset of the full expansion at any order. 



earized case, graphs do not contribute by their symmetry factors any more, as "insertion 
places" are missing. More precisely, in the so-called ladder-rainbow approximation we 
have r4(zi, ^2, -23, y) = ^{y, zi)5{y, Z2)6{y, z^) where the distribution d{x,y) = -^^S{x — y) 
is defined via 



dVy5{x,y)f{y) = f{x). 
Moreover, in this approximation the self-energy is given by 

nx,y) = \{-i\fG{x,y)^\x,y). 

As stated above, we also consider only the linearized version of the Dyson equation by 
replacing the exact propagator by the free one in the right hand side of (2.1). This means 
in particular that we do not sum over all self-energy chains. Finally we obtain the closed 
linear integral equation for the two-point function 



G{x,x') = A(x,x') 



6 



dVy^dVy^ A(x, yi)G{yi,y2)A^{yi,y2)A{y2,x') 



(2.3) 



which can be transformed into the differential equation [20]: 

\2 

(D^ + m2)(n^/ + m'^)G(x, x') A'^(x, x')G(x, x') = -i(D^ + m^)6(x, x'). 

6 

For maximally symmetric spacetimes and vacuum states this reduces to: 

A2 



{a + my -'—A\x,x') 



G{x, x) = — i(n -|- m )5{x, x). 



(2.4) 



In terms of the z variable and the hypergeometric operator H, using the formulas (1.1), 
the homogeneous ^ part of the preceding equation becomes: 



eA2]G(z) = 0, e = A76. 



(2.5) 



Note that for the non-linear Dyson-Schwinger equation, the corresponding differential 
equation would become highly non-linear even in flat space -see [15]-, and on top of that 



*The inhomogeneity proportional to a delta function will be responsible to a shift that can be absorbed 
in the integration constants 
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non-linearity we would have all the complexities of the hyper geometric operator on which 
we now focus. All further discussions of the non-linear case are postponed to future work. 

This ordinary differential equation is the central equation of our work and we will 
refer to it as the master equation in what follows. We will use the following boundary 
conditions, which are the natural generalization of the Bunch-Davies vacuum state (a 
direct derivation of these conditions from the integral equation will be given elsewhere): 
regularity for antipodal points (i.e at z = 0) and a flat spacetime singularity at short 
distances (at z — )• 1). The flat space singularity is proportional to (see Appendix B): 



^-2y^^+\ (2.6) 



Developing the relation z = cos^ {2k) ^^^^ /i — )• 0, we have the asymptotic behavior of G 
at z = 1: 



1 / Ve + 647r4 
' ' 4 V vr' 
where a is a given constant independent of e. 



G~^a(l-z)^ u = --\l^^^ +; 



3 Solution of the master equation 

We study the homogeneous master equation (2.5), which we reproduce here explicitly for 
convenience: 

^ 1 



This equation possesses one regular singular point at 2; = and two irregular singular 
points at z = 1 and z = 00. The local study performed around these points in Appendix 
C shows that the interaction term (proportional to e) is negligible near z = and is very 
important near z = 00. 

We then use a WKB approximation (or Liouville-Green for mathematicians) to obtain 
global approximations to the four fundamental solutions of (C.l). The solution that posses 
the correct z — t- 1 has the following IR behavior: 




^lni(-z) ^ 9TTy^- 



As the other solutions have a sub-leading contribution near z = 1, they (notably the 
one with u = 1 in the appendix C) can spoil the IR behavior by making it grow instead of 
decaying. This will be checked via a rigorous analysis elsewhere. 

We note that, while the approximation near z = 1 obtained via a certain WKB ap- 
proximation is only valid for strong couplings e ^ 1, the asymptotic behavior obtained via 
local analysis (see Appendix C ) is exact for any e. 

Finally an important remark: the phenomenon described here belongs to a large class 
of phenomena in mathematics and physics known as secular perturbation theory. The 
simplest such example is given by the boundary value problem [16]: 

y" + ey' + y = 0, for t > 0, y{0) = 0, y'{0) = 1. 
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To the first order in perturbation theory in e the solution is unbounded for large t: 

y{t) = sm{t) - ^etsin(i). 
Unlike this perturbative solution, the exact solution is bounded: 

y{t) = -^=l==e~-/^ sin (tVT^). 

In the case under consideration, the multiple scales method [16, 17] (consisting of intro- 
ducing a new scale t = et, treat it as independent of t, and use this freedom to kill the 
secular terms) yields, to the first iteration, the uniform approximation 

y(t) ~e-^*/2 gij^(^)^ 

We will investigate the applicability of these multiple scales methods to the study of strong 
IR effects in de Sitter in a future publication. Such secular behavior is usually related to 
a resonance in the system. Accordingly, we suggest that the expansion of the de Sitter 
background plays the role of a resonant external coupling. 

4 Dynamical mass generation? 

The authors of [7] retrieve the dynamically generated mass predicted by Starobinsky in 
[5]. Their approach is based on the constant "free" two-point function obtained in [4]. We 
show in this section how this result can be readily retrieved through our master equation. 

If we consider the constant two-point function Aq oc given in [4] (note that it is 
not really a free two-point function as it diverges for A — t- 0), our master equation (2.5) 
summing the ladder-rainbow diagrams has the four IR behaviors at z — t- oo: 

|^|(-v'9-4AAo-3)^ ^|(V9-4AAo-3) ^ ^|(-v'9+4AAo-3) ^ ^|(V9+4AAo-3) | _ 

The boundary conditions select ^2 (v^^~^^^~3) . As a mass m field two-point function (1.3) 
decays like z^i^^^^^^'^) ^ this implies a dynamically generated mass = XAq oc \/A. 



5 Discussion 



We end our work with a concise discussion of several important issues: 

• A crucial technical point in our analysis is that we are able to take full- advantage of 
de Sitter-symmetry by transforming the Dyson-Schwinger integral equation into an 
ordinary differential equation (depending only on the scalar variable z). 

• It is important to stress that our conclusions were drawn in a very restrictive context, 
namely the ladder-rainbow approximation, which is a drastic simplification of the full 
theory. A better control of this approximation has to be achieved, mainly by summing 
over a larger family of diagrams (for instance using the non-linear Dyson-Schwinger 
equation discussed before). We note however that the difficulty to rigorously justify 
the used approximation is a general trend in hard problems such as non-perturbative 
QFT, a classical example being the study of bound-states using the Bethe-Salpeter 
equation [21]. 
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• We have been able to avoid any discussion about UV-renormalization because we 
transform the integration over loops into differential equations. Hence renormaliza- 
tion only intervenes through the integration constants. This is actually the generic 
situation for the study of solutions of Dyson-Schwinger equations: for a kinetic renor- 
malization scheme, renormalization condition amount to fixing the boundary condi- 
tions of the equations. 

• The secularity of perturbation theory exists already in flat space. However, unlike the 
de Sitter case, it does not change the IR behavior of the two-point function, because of 
a rapidly decaying overall factor. Indeed the ladder-rainbow two-point function 
in flat space reads (see appendix B): 



• The Fourier transform and power-spectra of the non-perturbative functions G has to 
be studied. 

• The IR behavior of the graviton field in de Sitter spacetime is one of the most im- 
portant open questions in cosmology. We believe that the non-perturbative effects 
exhibited in the present work for the mmc scalar field can illuminate this issue. On 
top of these non-perturbative effects, we also expect an interesting interplay with 
non-trivial gauge artifacts similar to the ones already exhibited for the photon field 
in de Sitter in [18]. 

• In our study we neglect the back-reaction of the quantum fields on the background 
metric. However the fact that the interacting two-point function G decays in the IR is 
of paramount relevance to the difficult but crucial issue of whether the back-reaction 
can be neglected. Ultimately, this is equivalent to understanding the stability of the 
quantum de Sitter spacetime. 
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A Conventions 



We use the conventions of [19]. The metric is the mostly minus one. Consider the action 



S = d^xyj—g 



where g is the determinant of the metric. The free Feynman propagator {T(f)[x)(f)[x')) 
/\{x,x') verifies: 

{U^ + m^)A(x, x') = -i ^^^Z^\ 
The conventions for the Fourier transform are 
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A(x) 



ikx 



Feynman rules are given by: 



Propagator: A{x,x') 

Vertex: — iXs 

Mass renormalization: — i{mB — m'^ 

Field renormalization: — i{Z — l)n\x 



B Flat space master equation 

In the flat space the master equation 
solved and gives the four fundamental solutions with 



In the flat space the master equation (2.5), with A = 1/^ and A.{fi) = ^^^-i is readily 



/ Ve + 647r4 

The boundary conditions, included in the integral equation for instance, select the solution 
for which p = -^2^^^^. 

C Local analysis 

We study the homogeneous master equation (2.5), which we reproduce here explicitly for 
convenience: 

This equation possesses one regular singular point at 2 = and two irregular singular 
points at z = 1 and z = 00. We perform here a local analysis near these points. 
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C.l Frobenius series at z = 

The indical polynomial is Po{\) = A^(A^ — 1). Our equation being fourth-order and the 
roots of Pq being separated by integers, some care has to be taken in the application of the 
Frobenius method (see [22] for the explicit algorithm). We obtain the local behaviors: 

— , In z, 1, z 

z 

More precisely the first terms of the four solutions are given by: 
2z^ ( e 5' 



Inz ~ -, 

z 



G.{^) - ^ + 84 - (2I6 + ^) . + [-48 + (36 + ^) 



G^{z)^9z+(— ^ + 

^ ' V 9 46087r4 ' 



Inz ~ 21nz. 



This analysis implies that, the boundary condition imposing regularity at 2; = eliminates 
2 of the 4 solutions. 

C.2 Asymptotic behavior near z = 1 

e = 0. In this case the point z = 1 is a regular singular point. The free master equation 
being invariant under the transformation z — )• 1 — z, the asymptotic behaviors at z = 1 
can be directly obtained from the behavior near z = (also the equation is then exactly 
solvable) and are given by: 

l-z,-^,ln(l-z),ll. 



e > 0. Even if the point z = 1 is an irregular singular point, this "irregularity" appears 
only at the first order when z — t- 1. At zeroth-order in 1 — z, one can thus still define an 
indicial polynomial which is given by: 



2567r4 ■ 
The roots of which read: 



, ^1 / Ve + 647r4 

' = ^4V^^^' 
and are not separated by integers for generic e. The leading behaviors near 2; = 1 is thus 
given by 



n ^A X ^1 / Ve + 647r4 
We can verify that this is the full leading behavior at z = 1. 
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C.3 Asymptotic behavior near z ^ oo 

By making the ansatz G = and using the dominant balance method recursively [17], we 
are able to derive the leading behavior near infinity: 

LO^elJ{-z) 9TTy/\n{-z)\] ^_ 

1= 1 1= 5 ^ J-' 

3^/27r V2uj^ J\ 

Obtaining this result, necessitates some lengthy calculations that we choose not to repro- 
duce here. 



Guj 3 exp 

^3/2 J]^4 
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D WKB approximation 



Consider the general fourth-order differential equation (of course the parameter 5 can be 
introduced in many different ways, leading to more or less accurate approximations): 



54 



aoG = 0. 



When (5—7-0, each of the four elementary solutions in the WKB approximation read (in 
our situation, 5 = e"^/^, meaning that we perform a strong coupling expansion): 



G ~ exp 



1 



{So + 6Si + 6^82) 



where 



S()=uj 



00^ = 1, 



Si = ^[ 31ogao + 2 



^2 



03 



128a; 



/ao 



45a 



/ 2 



^ + 48a', + 



40a'(j 
ao 



+ 12a3^ - 32a2 



These solutions have the asymptotic behaviors 



(D.l) 



U\/e TT 



G, 



Ld Z— >00 



^3/2 ln4 (— z) 



exp 



w^ln2(-z) ^ 97rv/ln(-z) 



3^/27r 



V2u 



Comparing the WKB results to the exact local analysis performed before, we find that 
they are exactly the same as 2; — )■ 00. They also coincide near z — )• 1 in the strong coupling 
limit e 3> 1. This is a verification of our previous calculations. The "minimal" choice 
corresponding to the boundary value at z = 1 will be a; = —1, implying that the real part 
of the resummed two-point function decay at z — t- 00. 

However this WKB approximation breaks down near z = 0, preventing us from imple- 
menting the boundary conditions at z = 0. This breakdown can be explicitly seen because 
the following relations in the expansion G ~ exp (5*0 + 6 Si + 6'^S2)] do not hold near 
z = 0: 

5*0 > > S2. 

This was also expected, since the point z = is a degenerate regular singular point, i.e 
having equal exponents (see [23] p. 202). The global analysis of such equations, specially 
higher-order ones, is a non-trivial task (see [24] and references therein) that we will examine 
in a future work. 
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